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Abstract 

Many quasars and active galactic nuclei (AGN) appear in radio, op- 
tical, and X-ray maps, as a bright nuclear sources from which emerge 
single or double long, thin jets (Thomson et al., 1993). When observed 
with high angular resolution these jets show structure with bright knots 
separated by relatively dark regions. High percentages of polarization, 
sometimes more then 50% in some objects, indicates the nonthermal na- 
ture of the radiation which is well explained as the synchrotron radiation 
of the relativistic electrons in an ordered magnetic field. 

A strong collimation of jets is most probably connected with ordered 
magnetic fields. The mechanism of magnetic collimation, first suggested 
by Bisnovatyi-Kogan et al. (1969), was based on the initial charge sep- 
aration, leading to creation of oscillating electrical current, which pro- 
duces azimuthal magnetic field, preventing jet expansion and disappear- 
ance. Here we consider magnetic collimation, connected with torsional 
oscillations of a cylinder with elongated magnetic field. Instead of ini- 
tial blobs with charge separation, we consider a cylinder with a periodi- 
cally distributed initial rotation around the cylinder axis. The stabilizing 
azimuthal magnetic field is created here by torsional oscillations, where 
charge separation is not necessary. Approximate simplified model is de- 
veloped. Ordinary differential equation is derived, and solved numerically, 
what gives a possibility to estimate quantitatively the range of parameters 
where jets may be stabilized by torsional oscillations. 
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1 Introduction 

Objects of different scale and nature in the universe: from young and very old 
stars to galactic nuclei show existence of collimated outbursts - jets. Geometrical 
sizes of jets lay between parsecs and megaparsecs. The origin of jets is not well 
understood and only several qualitative mechanisms are proposed which are 
not justified by calculations. Theory of jets must give answers to three main 
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questions: how jets are formed? how are they stabilized? how do they radiate? 
The last question is related to the problem of the origin if relativistic particles 
in outbursts from AGN, where synchrotron emission is observed. Relativistic 
particles, ejected from the central machine rapidly loose their energy so the 
problem arises of particle acceleration inside the jet, see review of Bisnovatyi- 
Kogan (1993). 

It is convenient sometimes to investigate jets in a simple model of infinitely 
long circular cylinder, Chandrasekhar & Fermi (1953). The magnetic field in 
the collimated jets determines its direction, and the axial current may stabi- 
lize the jet's elongated form at large distances from the source (e.g. in AGNs), 
(Bisnovatyi-Kogan et al., 1969), see also Istomin and Pariev (1996), Beskin 
(2005). When observed with high angular resolution these jets show a struc- 
ture with bright knots separated by relatively dark regions (Thomson et al., 
1993). High percentages of polarization, sometimes exceeding 50%, indicate 
the nonthermal nature of the radiation, which is well explained as synchrotron 
emission of the relativistic electrons in an ordered magnetic field. Estimates 
of the lifetime of these electrons, based on the observed luminosities and spec- 
tra, often give values much less than the kinematic ages tk = d/c, where d is 
the distance of the emitting point from the central source. There is a neces- 
sity of continuous re-acceleration of the electrons in the jets in order to explain 
the observations. The acceleration mechanism for electrons in extragalactic 
jets proposed by Bisnovatyi-Kogan & Lovelace (1995), considers that intense 
long-wavelength electromagnetic oscillations accompany a relativistic jet, and 
the electromagnetic wave amplitudes envisioned are sufficient to give in situ 
acceleration of electrons to the very high energies observed > 10 13 eV. It was 
assumed that jets are formed by a sequence of outbursts from the nucleus with 
considerable charge separation at the moment of the outburst. The direction of 
motion of the outbursts is determined by the large-scale magnetic field. When 
the emitted wave is strong enough it washes out the medium around and the 
density can become very small, consisting only of the accelerated particles. The 
action of the oscillating knot is similar to the action of a pulsar, considered as an 
inclined magnetic rotator. Both emit strong electromagnetic waves, which could 
effectively accelerate particles. The model of enhanced oscillations of the cylin- 
der, and electromagnetic field around it, have been studied by Bisnovatyi-Kogan 
(2004). 

Here we consider stabilization of a jet by pure magnetohydrodynamic mech- 
anism, connected with torsional oscillations. Such type of oscillations in neu- 
tron stars had been considered by Bastrukov et al. (2002). We suggest that 
the matter in the jet is rotating, and different parts of the jet rotate in differ- 
ent directions. Such distribution of the rotational velocity produces azimuthal 
magnetic field, which prevents a disruption of the jet. The jet remains to be 
in a dynamical equilibrium, when it is representing a periodical, or quasiperi- 
odical structure along the axis, and its radius is oscillating with time all along 
the axis. The space and time period of oscillations depend on the conditions 
at jet formation: the length scale, the amplitude of the rotational velocity, and 
the strength of the magnetic field. The time period of oscillations should be 
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obtained during construction of the dynamical model, what also should show at 
which input parameters may exist a long jet, stabilized by torsional oscillations. 

2D nonstationary MHD calculations are needed, to solve the problem nu- 
merically. Here we construct a very simplified model of this phenomena, which, 
nonetheless, permits to confirm the reality of such stabilization, to estimate the 
range of parameters at which is takes place, and the connection between the 
time and space scales, magnetic field strength, and the amplitude of rotational 
velocity. 



2 Axially symmetric MHD equations 

Axially symmetric MHD equations at J^, for the perfect gas with an infinite 
conductivity, are written, in cylindric coordinates (r, 0, z), in the form (Landau 
and Lifshits, 1982; Bisnovatyi-Kogan, 2001) 
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P = P(p,T), E = E(p,T), f = f(p,T), (14) 

Here v = (v r , Dtp, v z ) is the velocity vector, B = (B r , B v , B z ) is the mag- 
netic field vector, j = (j r , j v , j z ) is the vector of the electrical current, p, P, E, 
are density, pressure, and internal energy, respectively, / is the cooling function, 
due to photons or neutrino. 

We consider a long cylinder with a magnetic field directed along its axis. 
This cylinder will expand unlimitly under the action of pressure and magnetic 
forces, so no confinement will be reached. The limitation of the radius of this 
cylinder could be possible in dynamic state, when the whole cylinder under- 
goes magneto-torsional oscillations. Such oscillations produce toroidal field, 
which prevent a radial expansion. There is a competition between the induced 
toroidal field, compressing the cylinder in radial direction, and gas pressure, 
together with the field along the cylinder axis (poloidal) , tending to increase its 
radius. During magneto-torsional oscillations there are phases, when either com- 
pression or expansion forces prevail, and, depending on the input parameters, 
we may expect tree kinds of a behavior of such cylinder. 

1. The oscillation amplitude is low, so the cylinder suffers unlimited expan- 
sion (no confinement) 

2. The oscillation amplitude is too high, so the pinch action of the toroidal 
field destroys the cylinder, and leads to formation of separated blobs. 

3. The oscillation amplitude is moderate, so the cylinder survives for an 
unlimited time, and its parameters (radius, density, magnetic field etc.) change 
periodically, or quasi-periodically in time. 

Solution of MHD equations (fl"|)- (fT4")) could give, in principle, the answer 
about a correctness of the above scenario. It is reasonable nevertheless to try to 
find a simple approximate way for obtaining a qualitative answer, and to make 
a rough estimation of parameters leading to different regimes. 



3 Profiling in axially symmetric MHD equations 

We'll try to simplify the system of equations in such a way, that the resulting 
system should contain the most important property of the dynamical competi- 
tion between different forces, to check the possibility of the dynamical confine- 
ment .We use for this purpose a profiling procedure. Let us neglect the gravity in 
the the direction of the cylinder axis (z), and approximate density by a function 
p(t,z), suggesting a uniform density along the radius, cylinder density. The 
components of the velocity and magnetic field are approximated as 
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v r =ra(t,z), v v = rQ(t,z), v z =0; (15) 

B r = rh r {t,z), B v = rh v (t,z),B z = B z (t,z). (16) 
In this case the current components (|5j)- (fTTj) are written as 

cr dh v . cr dh r . ch v 
3r = ~^~dz~' 3lp =4^~dz~' Jz = ^- 
After neglecting velocity v z along the axis, we should omit the corresponding 
Euler equation ([3]), and the radial pressure gradient is approximated by the 
linear function 

dP P 

where the constant A ~ 1 is connected with the equation of state, P(t, z) is the 
pressure, R(t, z) is the radius of the cylinder. In the subsequent consideration 
we consider an adiabatic case, where the polytropic equation of state P = K p 1 
is considered instead of the energy equation (fT3")) . Neglecting the z derivatives 
in the Poisson equation (fT2")l . we obtain tp G = nGpr 2 . Substituting (fT5"]) - ([T5j) 
into the original system of the equations, we obtain for the profiling functions 
the following equations 
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It follows from (f2"3|) -(f25" f relations, representing conservation of mass, and mag- 
netic flux equivalent to freezing condition 



P R 2 = C m (z), B z R 2 = C b (z), B z = -^rU (26) 
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In our subsequent consideration the arbitrary functions will be taken as con- 
stants: C m (z) — C m , Cb{z) — Cb- The algebraic relations (|26|) may be used 
instead of any two equations from (f2"3"j) - ([2l))) . 



4 Equilibrium configuration and linear oscilla- 
tions 

To check the properties of the approximate system (fT9|) -(f25 |) . we consider linear 
oscillations of the equilibrium, infinite, self-gravitating cylinder with uniform 
magnetic field and rotation along its axis: ^jf 1 — 0, ^j^- = 0. In equilibrium 
state (index " 0" ) we have 



p 
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Linearizing equations (fT9")) - (f2"2"|) around the equilibrium state ([2~T1) . we obtain 
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Linearizing of (f23|) - (f25|) . using the polytropic equation of state P = K p 7 , and 
looking for a solution in the form ~ expi(fcz — cut) we obtain 
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Here small perturbation values are taken without " 0" . From the last two equa- 
tions (f2"8")l we obtain using 
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First two equations (|28)l give after using (|29j) . ([30 
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Here we have introduced unperturbed sound (V s q), and Alfven (Vao) velocities, 
determined as 
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Equations (|31|1 lead to following dispersion equation 
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Let us consider several particular cases. 

1. Non-rotating, non-magnetized cylinder, Vao = = 0. It follows from 

CU) 
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Equation (|34p describes the only possible mode, remaining after fixing v z = 0, 
and uniform density over the radius. It corresponds to homologous mode of 
perturbation, where only radius and density (pressure) are oscillating. The 
frequency (|34j) becomes zero at 7 = 1 (isotherm). This degeneration is con- 
nected with the property of self-gravitating cylinder, which equilibrium is neu- 
tral (takes place at any radius) for the isothermal equation of state. This de- 
generation is equivalent to the well-known case at 7 = 4/3 for a spherical star 
(Chandrasekhar, 1939). While we are not interested in study of homologous 
oscillations, we'll take 7 = 1 in all father consideration. We have than from (|33[) 
the dispersion equation 
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which solution is written as 
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2. Rotating non-magnetized cylinder, Vao = 0, flo 7^ 0. Here the solution of 
dispersion equation uj 2 — $!§ ± £7§ describes the trivial mode cu 2 = 0, connected 
with pure rotational perturbations, and radial oscillations due to an action of 
the centrifugal force cu 2 = 2Hq (Fridman, Polyachenko, 1985). 

3. Non-rotating magnetized cylinder, Vao 7^ 0, ilo = 0. The solution of 
the dispersion equation uj 2 = k 2 V^ describes two different types of waves, 
propagating with Alfven velocity. To find out the nature of these oscillations 
note that (|3Tj) at 7 = 1, Sip = reduces to 
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The first type of wave corresponds to perturbation only of the radial velocity, 
a 7^ 0, £1 = 0. It follows that from ([3]), that h r ^ 0, h v = 0. This solution 
corresponds to the Alfven wave along the axis where only radius is perturbed, 
and no rotation appears. The second type of wave describes perturbations of 
angular velocity at constant radius, a = 0, fl ^ 0, where h r = 0, h v ^ 0. This 
wave corresponds to a pure torsional alfven wave along the cylinder axis. To 
obtain a standing wave ( e s ) we need a combination of two waves running in 
opposite directions (e r ±). 

a[e r+ cos (kz — Lot) + e r _ cos (kr + cut)] = 2a e s cos kr cos cot, 
h r [e r+ cos (kz — Lot) + e r _ cos (kr + tot)] 

= — B z q — [e r+ cos (kz — tot) — e r _ cos (kr + cot)] (38) 

w 

= —B z q — 2e s sinfcr sinwt = 2h r e s sinfcr sin cot. 

LU 

It is clear that radial velocity a and radial magnetic field h r are oscillating in 
opposite phases. For the torsional wave we have similar relations: 

O [e r + cos (kz — Lot) + e r _ cos (kr + cot)] = 2il e s cos kr cos cot, 
h v [e r+ cos (kz — uot) + e. r _ cos (kr + cot)] 
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CO 
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LU 

The rotational velocity and azimuthal magnetic field are also oscillating in op- 
posite phases. In reality we may have a mixture of these two degenerate modes 
where all perturbations are not zero, what is always takes place in nonlinear 
case. We see here that the approximate system of equations describes correctly 
small perturbations, connected with radial and torsional modes, so we expect 
that these modes will be described correctly in general nonlinear case. 



5 Farther simplification: reducing the problem 
to ordinary differential equation 

While in the relativistic jet the self-gravitating force is expected to be much 
less than the magnetic and pressure forces, we neglect gravity in the subsequent 
consideration. Without gravity the equilibrium static state of the cylinder does 
not exist. We need to solve numerically the system of nonlinear equations (fT§)) - 
(|2"2"|) . (f2"5]) , (|2T))) to check the possibility of the existence of a cylinder, which 
radius remains to be finite due to torsional oscillations (dynamic confinement). 

We'll try instead to reduce the system to ordinary equations, making addi- 
tional simplifications. Let us consider axially symmetric jet moving along z-axis 
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with a constant bulk motion velocity, in the comoving coordinate frame. In the 
jet which confinement is reached due to standing magneto-torsional oscillations, 
there are points along z-axis where rotational velocity always remains zero in 
this frame. Let us take f2 = in the plane z — 0. Let us consider standing wave 
torsional oscillations with the space period along z axis equal to zq. Then nodes 
with SI = are situated at z = ±n^, n = 0, 1, 2, .... Let us write the equations, 
describing the cylinder behavior in the plane z = 0, where f2 = 0. All values in 
this plane we denote by (~). We take also for simplicity A = 1. We have than 
equations in the plane z = as 



da i 2 _ K_ i C fc (dh r \ h%R 2 



dt B? 47rC ro \dz J z=0 2-kC, 
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from the equation (|2"2"]) ; 

^ = aA (44) 
dt v 1 

from the equation ([25)1 . The integrals of motion ([26]) in the plane z — are 
written as 

pi? 2 = C m , B z R 2 = C b , B z = ^p. (45) 
Initial conditions for the system (|4"0"[) - ([43|1 are 

^ = i?o, p = po = -S^i Bz = a = h r = h 9 = at t = 0. (46) 
R Q R 

In ([4*01) - (|43|) we have used relations 

'-"i^-^i- (47) 

valid for any time. If the cylinder rotational velocity is antisymmetric relative to 
the plane z = 0, fi = 0, and cylinder density distribution is symmetric relative to 
this plane, then we have extremum (maximum) of the azimuthal magnetic field 
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htj), with (^-iff-^j — 0, and zero value of h r = 0, which reaches an extremum 

(minimum) in this plane with (^p) z _ = 0. The product ap also reaches an 

extremum in the plane z = 0, so that ( ®12lLELA ] — o. The term with z 

. . \ dz Jz=0 

derivative in the equation (I43p is not equal to zero, and changes periodically 

during the torsional oscillations. We substitute approximately the derivative 

d/dz by the ratio 1/zo, where zq is the space period of the torsional oscillations 

along z axis. While Q — in the plane z = 0, its derivative along z is changing 

periodically with an amplitude f2o, and frequency w, which should be found 

from the solution of the problem. We approximate therefore 

= icos., (48) 
V dz J z=0 Za R2 

Finally, we have from (|4*0|) . (|4*5j) . (|4*4"|) the following approximate system of equa- 
tions, describing the non-linear torsional oscillations of the cylinder at given zq 
and Ho- 
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The combination of last two equations gives 
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dt Zq 

with the solution, satisfying initial condition (|46p . in the form 

- R 2 = CbS}o s . nujt (51) 

ZqUJ 

With account of ([5T|) the first and third equations in (j49|) give the equation 
~d(aR) /Cf,f2o \ 2 sin 2 tot 



V z w / 27rC m v ; 



Two differential equations, (|52|) and the third equation (|49|) determine the be- 
havior of the cylinder during magneto-torsional oscillations. Solutions where 
the radius does not go to infinity with time, determine a dynamically confined 
cylinder. Formation of blobs occurs, when the radius tends to zero. Long dy- 
namically confined cylinder exist when its radius is changing with time between 
two finite values. 
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6 Numerical solution 

Introduce non-dimensional variables 



R ' aR K y/K 

r = wt, V = tt> z = a ° ~ 77H2 ~ w ' ^0 = — — , (53) 

in which differential equations have a form 

dy dz 1 . „ r, „ , „ , , s 

-f- = z, — = - 1 - DsinM; y(0) = 1, z = at r = 0. (54) 
dr dr y 

Therefore, the problem is reduced to a system (|54[) with only two non-dimensional 
parameters D = (a*)', and ,(0), and the second one is taken equal 

to unity in farther consideration. All qualitatively different solutions are repro- 
duced inside this restricted set of parameters. Solution of this nonlinear system 
changes qualitatively with changing of the parameter D. 

The solution of this system was obtained numerically for D= 2, 2.1, 2.11, 
2.15, 2.2, 2.25, 2.28, 2.4, 2.5, 2.6, 2.9, 3.0. Roughly the solutions may be divided 
into 3 groups. 

1. At D < 2 there is no confinement, and radius grows to infinity after 
several low-amplitude oscillations (see Fig.l). 

2. With growing of D the amplitude of oscillations increase, and at D = 2.1 
radius is not growing to infinity, but is oscillating around some average value, 
forming rather complicated curve (Figs. 2-4). 

3. At D— 2.28 and larger the radius finally goes to zero with time, but with 
different behavior, depending on D. At D between 2.28 and 2.9 the dependence 
of the radius y with time may be very complicated, consisting of low-amplitude 
and large-amplitude oscillations, which finally lead to zero. The time at which 
radius becomes zero depends on D in rather peculiar way, and may happen 
at t < 100, like at D—2A, 2.6 (Figs. 7,10); or goes trough very large radius, 
and returned back to zero value at very large time r ~ 10 7 at D=2.5 (Fig. 
9). Starting from D = 3 and larger the solution becomes very simple, and 
radius goes to zero at r < 2.5 (Fig. 13), before the right side of the second 
equation (154j) returned to the positive value. The results of numerical solution 
are represented in Figs. 1-13. 



7 Discussion 

Let us consider, a jet which has an equation of state in the form P = Kp = v^p, 
v s — ° 2 /3> v s is the sound speed in the matter. For ultrarelativistic pair-plasma 
we have P = c 2 /3. The non-dimensional parameter D, as a function of the 
characteristic radius Rq, periodic length zq along z axis, initial density po, and 
magnetic field B z q, Qq and oj is written in the form 

D= j_^m_ (55) 
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The amplitude of oscillations fl, and u> should be found from the solution of 
the nonlinear system (fl9|) - ([25|) . together with determination of the interval of 
values of D at which confinement happens. In the approximate system (|54p only 
one parameter D characterizes different regimes, what for given values of Rq, 
zq, po, and B z q determines a function flo(D, uj) for the collimated jet. To find 
approximately a self-consistent model with ^q,uj(D) we may use the solution 
of linearized equations (|19|) - (f25|) with to ~ kV A from ([37]). The frequency of 
non-linear oscillations is smaller, and we may write 



2tt . , 

uj = a n kVA, a n < 1, k — — , (56) 



so that 



ttR 2 

.2 „2 (.2 1/2 „2 7lJD zO /r 7 \ 

= a n k V A = a n T . (57) 

Po z o 



Using it in (|55|) . with account of |53|) . we obtain 
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n 2 n2 o 2 n 2 2 , 2 D 2 ^ ,2 PO^s , r v . 

On the edge of the cylinder the rotational velocity cannot exceed the light 
velocity, so the solution with initial conditions in (|54[) . corresponding to yo = 1, 
has a physical sense only at < 2n i Da -i ~ j§^t- Taking a 2 = 0.1 for a 

strongly non- linear oscillations we obtain a very moderate restriction v s q < 2-. 
While in the intermediate collimation regime the outer tangential velocity is 
not changing significantly, this restriction would be enough also for the whole 
period of the time. To have the sound velocity not exceeding c/2, the jet should 
contain baryons, which density po cannot be very small, and its input in the 
total density in the jet should be larger than about 30%. 

The confinement by torsional oscillations starts at D — 2.1, and at D > 2.28 
the jet is divided into separate blobs according to Figs. 9-17. So, the confine- 
ment by magneto-torsional oscillations can be realized in the physically available 
situation, what is the main conclusion of this paper. 
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Figure 8: Time dependcnctkOf non-dimensional radius y (upper cOte), and 
non-dimensional velocity z (lower curve), for D = 2.5. 
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